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Abstract 

In this paper, we present some new results about the approximation of the Vlasov-Poisson 
system with a strong external magnetic field by the 2D finite Larmor radius model. The 
proofs within the present work are built by using two-scale convergence tools, and can be 
viewed as a new slant on previous works of Frenod & Sonnendriicker and Bostan on the 
2D finite Larmor Radius model. In a first part, we recall the physical and mathematical 
contexts. We also recall two main results from previous papers of Frenod & Sonnendriicker 
and Bostan. Then, we introduce a set of variables which are so-called canonical gyrokinetic 
coordinates, and we write the Vlasov equation in these new variables. Then, we establish 
some two-scale convergence and weak-* convergence results. 

1 Introduction 

Nowadays, domestic energy production by using magnetic confinement fusion (MCF) tech- 
niques is a huge technological and human challenge, as it is illustrated by the international 
scientific collaboration around ITER which is under construction in Cadarache (France). 
Since magnetic confinement, needed to reach nuclear fusion reaction, is a very complex phys- 
ical phenomenon, the mathematical models which are linked with this plasma physics subject 
need to be rigorously studied from theoretical and numerical points of view. Such a work 
programme based on rigorous mathematical studies and high precision numerical simulations 
can bring some additional informations about the behavior of the studied plasma before the 
launch of real experiments. 

The present paper can be viewed as a part of the recent work programme about the 
mathematical justification of the mathematical models which are used for numerical simula- 
tions of MCF experiments. Indeed, the first tokamak plasma models have been proposed by 
Littlejohn, Lee et al., Dubin et al. or Brizard et al. (see [5T], jTHJ [2D] , [E]> EHUD nevertheless 
most of these models were established by using formal assumptions. For ten years, many 
mathematicians have been working on mathematical justification of these models, especially 
the gyrokinetic approaches like guiding-center approximations and finite Larmor radius ap- 
proximations: many results in this research field are due to Frenod & Sonnendriicker et al. 
[TDJ HH E3 113 j Golse & Saint-Raymond [TH1IT5], Bostan [3] or, more recently, Han-Kwan 
[18]. These mathematical results mostly rely on two-scale convergence theory (see Allaire 
[3], Nguetseng [25]) or compactness arguments. 

In this paper, we are focused on the 2D finite Larmor radius model and its mathematical 
justification: more precisely, the goal is to make a synthesis of previous mathematical proofs 
of the convergence of (f e , E e ), where 

/ e (x, k, a, t) = f e (xi — V2k sin a, x 2 + v2fc cos a, V2k cos a, V2k sin a, i) , (1.1) 
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and where (f e , E e ) is the solution of the following 2D Vlasov-Poisson system 

d t fe + • Vx/ £ + (E £ + X - ( _^ ) ) • Vv/ £ = 0, 

< / e (x,v,0) = /°(x,v), (1.2) 

- Vi0 c = E £ , -Ai0 e = / f t dw -n e , 
\ Jr 2 

towards the couple (/, E) which is the solution of the 2D finite Larmor radius model given in 
[5] (see also Theorem [5] below) . The main results on this model are due to Sonnendrucker, 
Frenod and Bostan, and indicate that (/ e ,E e ) somehow weak-* converges to the solution 
of the finite Larmor radius model. However the proofs within these articles are based on 
various assumptions and use various tools. 

The first part of the present paper is devoted to a state-of-the-art about the two- 
dimensional finite Larmor radius approximation. Firstly, we recall the procedure which 
allows us to obtain the dimensionless model (jl.2l) from the complete Vlasov-Poisson model 
by considering specific assumptions. Then we recall the two-scale convergence theorem of 
Frenod & Sonnendrucker [T3] on the one hand, and the weak-* convergence theorem of 
Bostan [5 on the other hand. 

In a second part, we introduce a set of variables which are so-called canonical gyroki- 
netic coordinates and we reformulate the Vlasov-Poisson system (|1.2[) in these new vari- 
ables. Then, we establish a two-scale convergence theorem which only relies on Frenod & 
Sonnendriicker's assumptions. Finally, we deduce directly from this result Bostan's Finite 
Larmor radius model (see [5]) by adding a somehow non-physical assumption on the electric 
field which consists in considering a strong convergence of the sequence (E e ) e> o- 

2 State-of-the-art 

2.1 Scaling of the Vlasov-Poisson model 

This paragraph is devoted to the scaling of the following Vlasov-Poisson model: 
f t / + v-Vi/+— (E + vx B) ■ Vv/ = 0, 

J /(x,v,0) = / o (x,v), (2 .3) 
-Vx0 = E, -A x 0=— / fdv-— h e , 

where x = (2:1,2:2, $3) G K| is the position variable, v = (ii, v 2 ,v 3 ) £ R| is the velocity 
variable, t G R + is the time variable, / = /(x, v,£) is the ion distribution function, fi e is the 
electron density, E = E(x, t) is the self-consistent electric field generated by the ions and 
the electrons, B = B(x, t) is the magnetic field which is applied on the considered plasma, 
4> = <p(x,t) is the electric potential linked with E, e is the elementary charge and is the 
elementary mass of an ion. 

In this model, the external magnetic field B is assumed to be uniform and carried by the 
unit vector e3. We also assume that the electron density fi e is given for any (x, t) £ ]R| x M+. 
Following the same approach as in Bostan [S] , Frenod et al. [TU1 E] , Golse et al. [13 HE] 
and Han-Kwan [18] , we add the following assumptions: 

(i) The magnetic field is supposed to be strong, 
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(ii) The finite Larmor radius effects are taken into account, 

(iii) The ion gyroperiod is supposed to be small. 

We define the dimensionless variables and unknowns x' = (xi,x 2 ,x 3 ), v' = {v[,v 2 ,v 3 ), 
t', /', E' and & by 

Xi = L±x[, x 2 ~L±x' 2} X3 = L\\x' 3 , t = tt', \r = v\r', 
x'^vv^vv^vv^tt') =J f'ix^x^x'^v'-L^^v^t') , 

(2.4) 

E(Lj_ x[,L± x 2 ,L\\x 3 ,tt') = E'E'{x' l ,x' 2 ,x l 3 ,t') , 
4>(L±x[,L± x' 2 ,L\\x' 3 ,tt') = <j)4>'(x[,x 2 ,x' 3 ,t') . 

In these definitions, L± is the characteristic length in the direction perpendicular to the 
magnetic field, L|| is the characteristic length in the direction of the magnetic field, v is the 
characteristic velocity and t is the characteristic time. We also rescale the electron density 
as follows: 

n e {L^x x ,L^x 2 ,L\\ x 3 ,tt') = nn'^x^, x 2 , x 3 , t') . (2.5) 

Following the assumptions on the magnetic field B, we set B as being such that 

B = Ie 3 . (2.6) 

Then, we set L|| as the size of the physical device in game. We also link /, E and <fi with 
the characteristic Debye length Ad by 

— n — Xnen _ Ar> 2 en 

/ = =3> E =^— > ^~ > 2 " 7 

v e £o 

and we take Ad as the characteristic length in the direction perpendicular to the magnetic 
field, i.e. 

Ll = X^. (2.8) 

Since we want to take into account the smallness of the gyroperiod and the finite Larmor 
radius effects, we define the characteristic gyrofrequency TUl and the characteristic Larmor 
radius TZ as 

eB v 
uTi = , TZ=— . 2.9 

With these notations, the Vlasov-Poisson system is rescaled as follows: 



dff + tv 



Ad 1 \ ,7TT7„2 



'§v> 2 I . Vi ,/'+(^^E'+^v'xe 3 Vv,,/' = 



L \\ d 



/'(x',v',0) = /°'(x',v'), (2-10) 
E'= I , -A {s , iS ,J> - J^dl^' = / f'dv' 
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Taking into account the finite Larmor radius effects consists in considering a regime in 
which the Larmor radius is of the order of the Debye length. This implies 
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(2.11) 



Since the magnetic field is assumed to be strong, the Larmor radius is small when compared 
with the size of the physical domain. Then it is natural to take 



(2.12) 



where e > is small. 

Assumption (iii) can be translated in terms of characteristic scales by 



tuii 
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(2.13) 



Assumption (i) means that the magnetic force is much stronger than the electric force, so 

we consider 

Ee 

^^ = e. (2.14) 

v mi u>i 

Then, removing the primes and adding e in subscript, the rescaled Vlasov-Poisson model 
writes 
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dth + - . ,. 

e V v 2 



V (5l ,£ 2 )/ 6 + S 3 %,/e + (E e + x e 3 ) • Vj e = , 



/ e (x,V,0) = / e U (x,v) 



(2.15) 



-A 



f e dv - n e 



which is the model studied in previous works of Frenod & Sonnendrucker [14] . Golse & 
Saint-Raymond [HI [TC] , and Bostan [S] . 



2.2 Previous results 

In this paragraph, we recall two main results about the asymptotic behavior of the sequences 
{fe)e>o an d (E e ) e >o when e goes to 0. The first one is based on the use of two-scale con- 
vergence and homogenization techniques developed by Allaire [3] and Nguetseng [23], and 
was established by Frenod and Sonnendrucker in [14] . The second one relies on compactness 
arguments and was proved by Bostan in [5j. After recalling these two results, we discuss the 
main differences between them. These differences are the source of the motivation of the 
present paper. 

In order to simplify, we consider that the whole model (12.151) does not depend on X3 nor 
V3, and we assume that /f = f° for all e. Then, it is reduced to a singularly perturbed 2D 
Vlasov-Poisson model of the form 

dtfe + • Vx/ e + (E e + ~ ( _^ ) ) • Vv/ e =0, 

< / e (x,v,0) = /°(x,v), (2.16) 

- Vi0 £ = E £ , ~Ai0 e = / fedv-fie, 

{ Jr 2 
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where x = (xi,X2) & R 2 and v = (vi,vs) G 



The following theorem can be attributed to Frenod and Sonnendriicker |14j (see Theorem 
1.5) even if the setting of this paper is a charged particle beam Vlasov-Poisson model not 
involving any electron density. Nonetheless, the proof of [2] works again in the setting of 
model (|2.16l) which involves an electron density. 

Theorem 1 (Frenod & Sonnendriicker [E]). We assume that, for a fixed p > 2, f° is in 
L 1 (M 4 ) (1 L P (R 4 ), is positive everywhere and such that 



|v| 2 / (x,v)dxdv < +oo. (2.17) 
We also assume that h e does not depend on t, is in L 1 (R 2 ) l~l L 3 / 2 (R 2 ) and satisfies 

(2.18) 



/ (x, v) dx dv = / n e (x) dx 



Then, the sequence (/ £ ,E £ ) £>0 is bounded L°° (0, T; L p (R 4 )) x (L°°(0,T; VF 1:3 / 2 (R 2 ))) 2 in- 
dependently of e. Furthermore, by extracting some subsequences, 



-> F = F(x,v,T,t) two-scale in L°° (0, T; L|?(0, 2tt; Z7(R 4 ))) ; 

£ = £(5i,T,t) two-scale in (l m (0, T; L~(0, 2tt; W^ 3 / 2 (R 2 ))) 

Moreover, F is linked with G = G(y,u,*) € L°° (0, T; L p (R 4 )) fey i/ie re/ation 

F(x, v, r, t) = G(x + ft(-r) v, J?(-r) v, t) , 
and (G, £) is the solution of 



(2.19) 



(2.20) 



dtG(y,u,t) 



H(-a)8(y + K(a) u, a, t) da 



/ R(-a)£(y + 11(a) u,a,t) da 
Jo 



G(y,u,0) = ±-f°(y,u), 

Z7T 



£(x,T,t) = -Vi$(x,r,i). 



VyG(y,a,t) 
•V G G(y,u ) t)=0, 



(2.21) 



- A*$(x, r, i) = / G(x + ft(-r) v, JZ(-t) v, i) dv - — n e (x) : 



where 



K(t) 



sin r 1 — cos t 
cos t — 1 sin t 



JZ(r) 



cos r sm r 
— sin r cos r 



(2.22) 



In this theorem, L°£(0, 2tt; L p (R 4 )) stands for the space of functions h = h(x, v, r) being 
in L°°(0, 2tt; L p (R 4 )) and 27r-periodic with respect to r. 

As a consequence of this theorem, extracting some subsequences, we have 



U A / ini°°(0,T;L p (R 4 )), 

E £ A E in (i oo (0,T;W rl ' 3 / 2 (R 2 ))) 2 , 



(2.23) 
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where 



/(*,v,t) 



2;r 



F(x, v, r, i) <ir , and E(x, t) 



2tt 



£(x, r, t) dr . 



(2.24) 



By using the relation between F and G, we can easily remark that (/, E) is solution of 



d t f(x,v,t) 



2tv 



2tt 



2tt 



1Z(t — a) £ fx + 1Z((T — r) v, a, t) da 



2tt 



i?(r — cr) £ (x + 7?.(o- — t) v, cr, r, ) da 

/(x, v, 0) = i- y * /° (x + 7L{-t) v, fl(-r) v) dr , 

-V 4 ^(x,t) = E(x,t), -Ag0(x,t)= / /(x,v,i)dv-ne(x) 



VxF(x,v,r,t) dr 
VvF(x,v,r,t) dr = 0, 



(2.25) 

We notice that these equations still involve F and £. In order to make these dependencies 
disappear, Bostan has proposed in [S] a reformulation of the Vlasov equation in guiding- 
center coordinates. Before presenting it, we introduce the sequence (/ e ) e >o defined by 



/ e (x,v,i) = /Jx 



U2 
-Vi 



and, in the same spirit, we define the initial guiding-center distribution /° by 



/ u (x,v,i) = / u x + 



''2 
-«1 



(2.26) 



(2.27) 



Theorem 2 (Bostan [5]). FFe assume that fi e — 1, and £/ia£ / is 2it -periodic in x\ and £2, 
is positive everywhere and satisfies 



2ir /.27T 



Jo 



/°(x,v)dxdv = 1, 



2-7T /.2lT 



Jo 



v| 2 /°(x,v)dxdv < +00. (2.28) 



VFe a/so assume that there exists F° € L°°(R + ) (1 L 1 (IR.+ ; r dr) sttc/i i/ia£ 

V (x, v) e [0, 27r] 2 x R 2 , /°(x,v) < F°(|v|) . (2.29) 

We also assume that (E e ) e > admits a strong limit denoted with E in [L 2 (0, T; L^.([0, 27r] 2 ))) 2 . 
Then, up to a subsequence, f e weakly-* converges to a function f = /(x, v, t) in L°°([0,T) x 
R 2 ;L-([0,2tt] 2 )) verifying 

/(x,v,t) = ^- 5 (x,^- ! t), (2.30) 
where g — g(x, k, t) is the solution of 
' d t g + (£ 2 ) d Xl g - (£1) d X2 g = , 

/•2tT 

g(x, fc,0)= / /°(x, V2fc cos a, Vzfc sin a) da , 
Jo 

1 f 271 

(£)(-x,k,t) = — / E(x! - V2k sma,x 2 + V2k cos a,t) da, (2.31) 
2-7T J 

- Vs0(5,t) =E(x,t), 

Aj0(x,t) = — / / + \/2fc sina, X2 — %/2fc cos a, fc,i) dadk — 1 . 



2^ Jo 
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Iii this theorem, L 2 ^ ([0, 27r] 2 ) stands for the space of functions h — h(x) being in 
L 2 ([0,27r] 2 ) and 27r-periodic with respect to x\ and x 2 . 

This last result introduces a mathematical justification of the approximation of the 
Vlasov-Poisson model (|2.16[) by the finite Larmor radius model which is exactly (I2.30[) - (I2.31[) . 
However, in order to prove this convergence result, Bostan considered stronger assumptions 
on /° and n e than needed to get existence of the weak-* limit (/, E) from Theorem Q] the 
initial distribution f° is supposed to be 27r-periodic in x\ and i 2 and (E e ) e> o is supposed 
to admit a strong limit in some Banach space. 



3 Convergence result in canonical gyrokinetic coordi- 
nates 

This section is devoted to a two-scale convergence result for model (|2.16|) written in a new set 
of variables which are so-called canonical gyrokinetic coordinates. Then, by adding a non- 
physical hypothesis for the electric field E e , we obtain straightforwardly Bostan's weak-* 
convergence result. 

3.1 Reformulation of Vlasov equation 

Following the ideas of Littlejohn [21] , Lee [El [20] , and Brizard et al. [H [7] , we define the 
variables (xi,x 2 ,k,a) el 2 x K + x [0,27r] by linking them with [x\, x 2 , V\, v 2 ) € R 4 by 

x\ = x\ - v 2 , vi = \/2kcosa, . 

x 2 — x 2 + vi , v 2 = V2k sin a. 

This set of variables is so-called canonical gyrokinetic coordinates: indeed, if we define the 
characteristics X±, X 2 , K, A linked with X\, Xa, k, a. by 

Xi = Xx-Vz, V x = V2KcosA, 

X 2 = X 2 + Vi, V 2 = V2KsinA, [ ' ' 

where X\, X 2 , V\, V 2 are the characteristics associated with the Vlasov equation (|2.16l a). i.e. 
satisfying 



d t Xi{t) = -V x {t), d t V x {t) = E etl (X 1 (t),X 2 (t),t) + -V 2 (t). 
e e 

d t X 2 (t) = -V 2 (t), 8 t V 2 (t) = E^(X 1 (t),X 2 (t),t)--V 1 (t) : 
e e 



(3.3) 



we have 



(3.4) 



d t X x {t) = -d X2 H t {X 1 {t),X 2 {t),K{t\A{t),t), 

d t X 2 (t) = d Xl H e (X x (t),X 2 {t),K{t),A(t),t) , 

d t K(t) = daH^X^X^Kit^Ait)^), 

d t A(t) = -duH^X^X^K^.A^t), 

where the hamiltonian function H t is defined by 

k 

H e (xi,x 2 ,k,a,t) = - + (f> e (xx,x 2 ,k, a,t) , (3.5) 
e 

and </> e is linked with </> e by the relation 

(f> e (xi, x 2 , k, a, t) = 4> e (xi — V2k sin a, x 2 + V2k cos a, t) . (3-6) 
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Then it is straightforward to see that, in the gyrokinetic canonical coordinates, the Vlasov- 
Poisson system (|2.16[) has the following shape: 



d t fc + E ta d x J e - E e ,i d x J e + V2k (E e>1 cosa + E^ 2 sina) d k f e 

E e , 2 cosa- E e<1 since 1 

4/e Q / e = 0, 



'2k e 
/ e (x, fc, a, 0) = /° (x\ — V2k sin a, 22 + cos a, V2fc cos a, \/2fc sin a) , 
E e (x, fc, a, t) = E e (xi — V2k sin a, £2 + V^fc cos a, t) , 
E e (5,t) - -V*&(x,t), 

/ e (xi + V2k sin a, x 2 — V2k cos a, fc, a, t) & da — h e (x, t) , 



(3.7) 



Jo 

where x = (xi,X2) and x = (act, X2), and where / e = / e (x, k, a, i) and E e = E e (x, k, a, t) are 
linked with f e and E c by 

/ e (x, k,a,t) — f e (xi — \/2fc sina, 22 + V2fc cosa, v2fc cosa, v2fc sina, t) , 

(3-8) 

E e (x, k,a,t) = iEiJxi — v2k sina,X2 + \2k cosa,t) . 

3.2 Two-scale convergence 

We set the following notations 

n = M 2 x K + x [0, 2tt] , r = R 2 x R + , S = R + x [0, 2tt] , (3.9) 
and we consider the following Banach spaces, involving periodicity with respect to a: 

L p # (0,2tt;L p (T)) = {/ e L p {Vt) : f is periodic ma}, 
W^' p (0, 2ir; W 1 '^)) = {/ G W^ p (ft), /(., ., 0) = /(., ., 2tt) : / is periodic ma}, 
Wf p (0, 2^; W 2 ' P (T)) = {/ e W 2 ' p (fi), /(., ., 0) = /(., ., 2tt), 5 Q /(., ., 0) = d a f(., ., 2tt) : 

/ is periodic ma}, 

and we can state the following theorem. 

Theorem 3. We assume that, for a fixed p > 2, f° and n e satisfy the assumptions of Theo- 
rem^ Then sequences (f e )e>0 an d (E c ) c >o of system \3. 7[ ) are bounded independently of e 
in L°°(0,T;LP (0,2n; L p (T))) and (L°°(0,T; W 1,3 ^ 2 (R 2 ))) respectively. As a consequence, 
there exist F — F(x, k, a, r, t) and £ = £(x, r, t) such that, extracting some subsequences, 

ft — > F two-scale m (0,T; (0,2tt; L p # (0,2ir; Lr(T)))) , 
E e — > £ two-scale in (L°° (0, T; L~(0, 2tt; VK 1 ' 3 / 2 ^ 2 )))) 2 . 

Furthermore, there exist G = G[x,k,a,t) € L°°{0, T; 2^(0, 2tt; L p (T))) and£ = £(x, ft, a,r,i) £ 
(L x {Q,T\L^{Q,2Tr\W^ 3,2 {Q,2^W l '^ 2 {V))))) 2 such that 

F(x, k, a, r, t) = G(x, fc, a + r, t) , (3-11) 
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£(x, fc, a,r, t) = £(xi — \p2M, sina,X2 + V2k cosa,r, t) , (3-12) 

and verifying 

( d t G + (£ 2 ) 8 X1 G - (£,) d X2 G+ (J- fe ) d k G + (T a ) d a G = , 

G(x, fc, a, 0) = — /°(xi — V%k sin a, x 2 + V%k cos a, \/2k cos a, \/2k sin a) , 

-Vi*(x,T,i) = £(x,T,t), 

— Aj$(x, r, i) = / G(xi + \/2fc sin a, £2 — \/2fc cos a, fc, a + r, i) dfcda (3.13) 
Js 

--n e (x), 



7-fe(x, fc, a, r, t) = v 2fc (£i(x, fc, a, r, t) cos a + ^(x, fc, a, r, f) sin a) 

£2(x, fc, a, r, t) cosa — £ i(x, fc, a, r, i) sina 
y- Q (x, fc, a, t, i) = 



V2fc 

where the notation (•) stands for 



r 2ir 

(u) (x, fc, a, t) = I u(x, k, a — r,T,t) dr . (3-14) 



(i 



Proof of Theorem [3J Several parts of this proof are only sketched since they can be redundant 
with [T3]. However, more details can be found in Mouton [23] , 

Following the same way as in |14j . we prove that, under the assumptions of Theorem [1] 
we have 

l|AM)IL* (0,2^)) = ll/V(*) < V ^ > ( 3 - 15 ) 

and, dehning p € as 

p e (Sc,t) = / / e (xi + \/2fc sina,X2 — V^A cosa, k,a,t) dkda , (3.16) 
Js 

that the sequence (/5 c ) c >o is bounded in L°° (0, T; i 3 / 2 (M 2 )) independently of e. Then, we 
deduce that (0 e )e>o and (E e ) e>0 are bounded independently of e in L°° (0, T; W 2 ' 3 ^ 2 (R 2 )) 
and (£°°(0, T; W /1 ' 3 / 2 (R 2 ))) 2 respectively. As a consequence, there exist F = F(x, fc, a, r, £), 
$ = $(x, r, f) and £ = £(x, r, i) such that 

U — > F two-scale in L°°(0,r;i^(0,27r;L^(0,27r;LP(r)))) , 

4> e — > $ two-scale in L oo (0,T;i~(0,27r;M/ 2 ^/ 2 (E 2 ))) , (3.17) 

E e — ► £ two-scalein (0,T;£^(0, 2tt; M/^ 3 / 2 (R 2 )))) 2 . 

Considering a compact set K C T, we easily remark that the sequence (</> e ) e > defined by 
(EU) is bounded in L°° (0, T; W^ ,3/2 (0, 2tt; W 1 ' 3 / 2 ^ ))) and that all its second order deriva- 
tives except dl<j>z are bounded independently of e in L°° (0, T; L 3 J 2 (0, 2tt; L 3 / 2 (A))). The 

sequence (E e ) e>0 defined by (|Ob1 is bounded in (A°°(0, T; W^' 3//2 (0, 2n; Fy 1,3 / 2 (A)))) 2 
independently of e. As a consequence, we claim that there exist $ = <f>(x, fc, a, t, i) and 
£ = £(x, fc, a, r, t) such that 



£ — -> $ two-scale in L°°(0, T; ^(0, 2tt; ^' 3/2 (0, 2tt; T4^ 1 ' 3 / 2 (A)))) , 
E e — > £ two-scale in (L°°(0,T; L^(0,2ir; W^ 3/2 (0,2ir;W h3 / 2 {K))))y 



(3.18) 
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Furthermore, we remark that $ and £ are linked with $ and £ by the formula 
£ (x, fc, a, r, t) = £(.Ti — \Flk sina, X2 + v2fc cos a, r, t) , 
$(x, fc, a, t, i) = — v2fc sin a, X2 + V2k cos a, r, t) . 
Then the vector function A £ defined by 

/ E e , 2 \ 



(3.19) 



A f = 



/ 


~d X2 (/)e 


\ 






9 Xl (f>e 








d a <t>e 






V 




J 





\ 



E e<2 

-Ee.l 

2k [E ti i cos a + E ti2 sin a) 
E t 2 cos a — E e i sin a 



(3.20) 



2A- 



has its three first components which are bounded in L°° (0, T; (0, 2n; W 1 ' 3 ^ 2 ^))) , in- 

dependently of e and its fourth one in L°° (0, T; i# /2 (0, 2tt; L 3 / 2 (iT))) and admits a two-scale 
limit denoted A = A(x,k,a,r,t) in (L°° (0, T; i^(0, 2tt; ^ 3/2 (0, 2tt; l^ 1 ' 3 /^^ ))))) 3 x 

L°° (Q,T; L^(0,2tt; L 3 J 2 (0,2tt; L 3 / 2 (K)))) . The convergence of the three first components 
is the consequence of classical embedding of Sobolev spaces in LP spaces. Concerning the 
convergence of the fourth one we need to use that E e ^i and E tL ^, are bounded independently 
of e in L°°(0,T',W^ 3/2 (0,2w, W^ 3 / 2 (K))) and consequently in L°°(0, T\ L^(0, 2tt; L e (K))) 
and that (2fc) -1 / 2 is, for any fc max and any q < 2, in L q (0, fc max ). The vector function A is 
linked with <!> and £ as follows: 

£2 \ 

-£1 

^2k (£ 1 cos a + £2 sin a) 
£2 cos a — £\ sin a 



A 







\ 


















V 




J 





\ 



(3.21) 



1 



N -V/, 



(3.22) 



In order to establish the two-scale limit model, we cannot simply apply Theorem 1.3 of 
[14] : indeed, the formulation (|3. 71 a) of Vlasov equation does not fit with the assumptions 

1 

which are needed for applying this theorem since the differential operator / 1— > — d a f cannot 

e 

be written under the form 

fx 1 \ 

X2 
k 

V a ) 

where M is a constant square matrix satisfying Tr(M) = 0, and N 6 im(M). However, the 
approach which is considered in [14] can be adapted to the present case. 

Firstly, we prove that there exists a function G such that F(x, k, a, r, t) — G(x, k, a+r, t). 
To reach such a result, we consider a test function -0 = tp(x,k,a,T,t) on x [0, 27r] x [0, T] 
which is 27r-periodic in a and r. If we multiply (|3.7l a) by 0(x, k,a, -,£) and integrate over 
fix [0,T], we obtain 

l-T 



f e (x,k,a,t) d t i/j(x,k,a,-,t) + -d T ^p(x,k,a,-,t) 

L v e ' e v e ' 

+ A £ (x, k, a,t) ■ V/ e (x, k,a,t) - - d a ip(x, k,a 7 -,t) 



efac d/c da 



(3.23) 



f°(xi — v2fc sin a, a?2 + v2fc cos a, v2fc cos a, v2fc sin a) 
x ip(xi, X2, k, a, 0, 0) dxi dx2 dk da . 
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Multiplying (|3 . 23[) by e and letting e — > 0, we obtain the weak formulation of d T F—d a F = 0, 
which indicates that there exists a function G G L°°(0, T; Lg(0, 2ir; L P (T))) such that 

F(x, k, a, r, t) = G(x, fc, a + t, i) . (3.24) 

Secondly, we introduce the sequence (<? e ) e >o defined by 

g e (x,k,a,t) = f e (x,k,a- -,tj . (3.25) 

In the spirit of [T3] , we prove that g e strongly converges to 2-7T G in a given Banach space. For 
that, we notice that, up to a subsequence, g c two-scale converges to G in L°°(0, T; (0, 2tt; 
L p # (0, 2tt; LP(T)))) since we have 



(3.26) 



/ / g e (x,k,a,t) ip(x,k,a, dx dk dadt 
Jo Jn e 

= / / / £ (x, fc, a, t) ip(x, k, a -\ — , -,t) dx dkdadt 
Jo Jn 6 e 

p2tz pT p 

— > I / / F(x, k,a,r,t) ip(x, k, a + T,T,t) dx dk dadt dr 
Jo Jo Jn 

p*2tt pT p 

= J / / G(x, k, a,t) ip(x, k,a,r,t) dx dk da dt dr , 
Jo Jo Jo 

for any test function tp on f2 x [0, 27r] x [0, T] which is 2-7r-periodic in a and r. As a consequence, 
g e weakly-* converges to 2tt G in L°° (0, T; (0, 2tt; L p (r))) . 

Let us prove that this weak-* convergence is a strong convergence in a given Banach 
space. This is the aim of the following lemma: 

Lemma 4. For any compact subset K ofT, and up to a subsequence, g t strongly converges 
to 2ttG in L°° (0,T;(W^ 3/2 (0,2tt;Wo' 3/2 {K)))*) . 

In this Lemma (W^' 3/2 (0, 2tt; W h3/2 {K)))* stands for the dual of VF^ 3/2 (0, 2tt; W„ 1,3/2 (K)). 

Proof of Lemma^ For any compact subset K of T, (<? e ) e> o and (A e ) e> o are bounded 
independently of e in L°° (0, T; L\ (0, 2tt; L p (iQ)) and (i°°(0, T; !F#' 3/2 (0, 2tt; Fy 1 - 3 / 2 ^)))) 3 
x_L°°(0, T; L^ 2 (0, 27r; L 3 / 2 (i4T))) respectively. Then, remarking that g e is solution of 

dtg € (x, k, a, t) + A c (x, k, a , t) ■ Vg e (x, k,a,t) = , 

e (3.27) 

# £ (x, fc, a, 0) = f°(xi — \[2~k sin a, x% + \J~2~k cos a, \[2~k cos a, V2k sin a) , 

we use similar arguments as the ones given after equation (|3.20l) to deduce 

1. (A e ) e>0 is bounded independently of e in (L°°(0, T; L| /2 (0, 2tt; L 3 / 2 (iT)))) 4 , 

2. (g £ (x,k,a,t) A £ (x,k,a- f,*)) e>0 is bounded in (L°°(0, T; L^(0, 2tt; Z7(if)))) 4 inde- 
pendently of e with r defined by £ = i + | (r € ]1, | [), 

3- (<9 t g e ) £>0 is bounded in L°°(0,T; (VF^' r * (0, 2tt; W Q hr ' {K)))*) independently of e with 
1. 



r* r 
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If r < 3/2, the embedding W^ r * (0, 2tt; W 1,r * (K)) C W^ 3/2 {Q,2tt;Wq 3/2 (K)) is compact 
with density. Furthermore, Rellich-Kondrakov's theorem (see p]) gives the compact embed- 
ding L^.(0, 27r; L p (K)) C (W^' 3/2 (0,27r; W 1,3/2 (i4T)))* since p > 2. Then, we apply Aubin- 
Lions' lemma (see [H]) and we prove that the functional space hi defined by 

U = [u e L°°(0,T;L!L (0,2tt:L p {K))) : 

1 , (3-28) 

d t ueL°°(0,T;(W^ r (0,2n;W^ r (*)))*)}, 

is compactly embedded in L°°(0, T; (W^ 2 (0, 2tt; Wq ' 3 ^ 2 (K )))*). Since g £ € U for all e, we 
deduce that the weak-* convergence of g e to 2ttG in L°° (0,T; L p # (0,2ir; LP(K))) is a strong 
convergence in L°°(0,T; (W^ 3/2 (0, 2tt; W 3 ,2 {K)))*) . 

If r > 3/2, the compact embedding (W^ r " (0, 2tt; W 1 ' r '* (X)))* C (W^' 3/2 (0, 2tt; W X ' 3/2 (#)))* 
is gotten directly. If we introduce the functional space W defined by 

W = {ueL™(Q,T-L p JQ,2>ir-,L p {K))) : 

, , , (3-29) 

d t u € L°°(0,T; (^ 3/2 (0,27r;<' 3/2 (X)))*)} , 

we remark that the sequence (g e ) c >o is bounded in hi' independently of e. By using Aubin- 
Lions' lemma, we prove that W C L°°(0,T; (W^ 3/2 (0, 2tt; W 1,3/2 (lf)))*) is a compact em- 
bedding. Then, g e strongly converges to 2tt G in L°° (0, T; (W^' 3/2 (0, 2tt; W 1,3/2 (#)))*) . □ 



dx rffc da dt 

(3.30) 



Let us finish the proof of Theorem [3] by establishing a transport equation satisfied by G. 
Let us consider a test function ip = ipipc, fc, a, t) on VL. We denote its compact support in T 
by X and we assume that ip and its first order derivatives are 27r-periodic in the a direction. 
Then we have 

III ge(x,k,ot,t) dtip(x.,k,a,t)dx.dk dadt 

JO JO J K 

+ / .9e( x > k, a, t) A c (x, k, a , t) ■ Vip(x, k, a,t) 

Jo Jo Jk 6 

+ / / f°( x i ~ sin a, x% + v2k cos a, \Flk cos a, \Flk sin a) 
Jo Jk 

x i/'( x j a, 0) d~x.dk da = . 

Since the vector function Z £ (x, fc, a, t) = <7 e (x, fc, a, t) A £ (x, fc, a — |,t) is bounded indepen- 
dently of e in (L°°(0,T;£^(0,27r;L r (ir)))) 4 with r e]l,§[, there exists a function Z = 
Z(x,k,a,r,t) in (L°°(0, T; L# (0, 2tt; L#(0, 2tt; L r (K))))) 4 such that Z e two-scale converges 

to Z in (L°°(0, T; L3f(Q, 2tt; LU0, 2tt; L r (iT))))) 4 and weakly-* converges to / Z(; •, •, r, •) dr 

Jo 

in (L°°(0,T;L%(0,2Tr,L r (K)))) 4 . 

If we now restrict our choice of test functions ip to the one having their support satisfying 

K n (R 2 x {0}) =0, (3.31) 

then, A e (x,fc,a- \,t) is bounded in (L°° (0, T; W^' 3/2 (0, 2tt; W 1 ' 3 / 2 ^)))) 4 independently 
of e and admits a two-scale limit in space (i°°(0,T; L%(Q, 2tt; W^' 3/2 (0, 2tt; W^ 3 / 2 (if))))) 4 
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which is exactly A Since sequence (g £ ) £>0 converges to 2it G in space L°°(0, T; (Wi ,3/,2 (0, 2it; 
Wq' S (K)))*) strongly, we can deduce that 



Z = 2nGA, (3.32) 

almost everywhere on f2 x [0, T]\{(x, k,a,t),k = 0} and then almost everywhere on f2 x [0, T], 
and we obtain the weak formulation of 



d t G(x, k, a, t) 



/ A(x, k,a — r, t, t) dr 
Jo 



VG(x, fc, a, t) = 0, 

(3.33) 

G(x, k, a, 0) = -^—f°(xi — \/2fc sina,X2 + \/2fc cos a, Vzfc cos a, \/2fc sin a) , 

when e — > in (|3.30p . Finally, we extend this conclusion to every compact subsets K of T 
by remarking that the measure of M 2 x {0} is null. 

In order to obtain Poisson type equations p. 131 c) and (|3.13l d). we consider a test function 
f/> = ?/>(x, r, t) on R 2 x [0, 2tt] x [0, T] which is 27r-periodic in r, we multiply E c (x, t), Vx</> e (x, t) 
and A x ^ e (x, i) by ipi*, f , i) and we integrate in x and t. We obtain 



[ ( E e (x,t)ip(x,-,t)dxdt^ [ [ [ £(x,T,t)ip(k,T,t)dxdtdT , (3.34) 
Jo Jr 2 e Jo Jo Jr 2 

[ [ Vi^ £ (x,t)V(x,-,t)dxdt-> [[[ V*$>(x,T,t)ijj(x,T,t)dxdtdT , (3.35) 
Jo Jr 2 e Jo Jo Jk 2 



'o 

/ / A^ e (x,t)ip(x,-,t)dxdt^ Iff Az${x,T,t)tp(x,T,t)dxdtdT , (3.36) 
Jo Jr 2 e Jo Jo Jk 2 

when e converges to 0. Since F is the two-scale limit of (f e ) e >o, we also have 



/ / I / / £ (^! + sina, X2 — v^fc cosa, fc, a, t) dkda J ^(x, -, t) dxdt 
Jo Jr 2 \Js J e 

— >• / / / / + V2fc sina, £2 — V / 2fc cosa, fc, a, r, i) 

Jo Jo Jr 2 Js 

x ^(x, r, t) <ifc da dx dt dr , 



(3.37) 



when e — > 0. Then, gathering convergence results (|3.34p - (|3.37p . we obtain the weak formu- 
lation of (|3TT51 c) and (|3TT51 d). □ 

3.3 Weak-* convergence under non-physical hypothesis 

Firstly, we have a direct corollary of Theorem [3] 
Corollary 5. Up to some subsequences, 

• f e weakly-* converges to f E L°° (0,T; L^(0,2vr; i p (r))) urai/i / denned &y 

f(x,k,a,t)= / F(x,k,a,r,t) dr , (3.38) 
Jo 

• E e weakly-* converges to E G (L°° (0, T; 1F^ 3 / 2 (R 2 ))) 2 luftfc E denned fey 

E(x,t)=/ £(x,T,t)dr, (3.39) 
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4> e weakly-* converges to G L°° (0, T; 1F 2:3 / 2 (IR 2 )) to^/i defined by 



271 



x,t) = / <f>(x,T,t)dT. (3.40) 



From now, we add the following non-physical assumption: 

E e — J- E strongly. (3.41) 
Lemma 6. Under the hypotheses of Theorem^ and \3.4-l\j , we have 



V(x,fc,a,t) e fi x [0,T] , $(x,fc,a-r,T,t)dr^ = 0. (3.42) 

Proof of Lemma\^ We consider a compact subset K of T and a test function ?/> = ip(x, k, a, t) 
on fix [0, T) which support in (x, k) is included in K. Since (<j> e ) e > o is bounded independently 
of e is L oo (0,T;^ 3/2 (0,27r;l^ 1 ^ 3 / 2 (i ; i:))) and that all its second order derivatives, except 
d 2 k (j) tl are also bounded in L°° (0, T\ L^ 2 (0, 2n; L 3 / 2 (K))) , we have 



r f Z7r f t 

III ^xd a (j> t (x,k,a ,t)ip(x,k 7 a,t) dxdkdadt 

Jo Jo Jr e 

,T ,2tt , t 

= — I / / V x ^e(x, k, a, t) {d a "4>){x, k, a H — ,t)dxdkdadt 
Jo Jo Jr e 

[ I I ' E e (x,k,a,t)(d a ilj)(x,k,a + -,t)dxdkdadt (3-43) 
Jo Jo Jr e 



o Jo Jr 

x (d a tp) (x\ + V2k sin a, X2 — V2k cos a, k, a H — , t) rfx dfc da dt . 
Since E e converges to E strongly, we obtain 

lim ((( V x <9 Q (/> e (x, k, a ,t)if)(x,k,a,t)dxdkdadt 

e ^°io Jo ir e 

/ / E(x,i) 
o Jo Jr 

2tt 



X 

T ,.2 



o Jo Jr 

2tt 





/ {d a ip) (i l + V2k sin a, X2 — V2k cos a, fc, a + r, t) dr dx dk da dt 
Jo 

E(x,t) 

y <9 T ^/>(5i + \/2fc sin a, £2 — V2k cos a, fc, a + r, t)^ dr dxdk da dt 



(3.44) 



= 0. 



On another hand, we deduce from (|3.18l) that we also have 

T f 2 * f t 

/ / V x <9 Q (/> e (x, k, a , t) ip(x,k,a,t) dxdk dadt 

Jo Jr e , 

(3.45) 



J J J ^x^a ( J ^(x, fc,a — T,T,t) drj ip(x, k,a,t) dxdk da dt . 
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Then it is straightforward that 



2jv 



V x d Q / $(x,k,a- T,T,t)dr =0, (3.46) 



o 



in L°° (0,T; l| /2 (0,2tt; L 3 / 2 (K))) and that (x,k,a,t) i-> d a ^ $(x, fc, a-r, r, t) dr^j does 
not depend on x. 

Following the same approach, we use the fact that the sequences (<9f a 4> e )t > o and (d 2 (j) e ) e > o 
are bounded in L°°(0, T; L^ /2 (0, 2tt; L 3 / 2 (K))) independently of e to obtain 



in L^(0,T;L^ /2 (0,27r;L 3 / 2 (iO)). On the other hand, we have 



T r 2ir 



JO 



t 



d k Q (/) e (x, fc, a , t) ?/>(x, fc, a, i) efocdfc da eft 



e 



a 27r f E e2 {x,k,a,t) cos a ^ w i 
/ — : = (d a ip)(x.,k,a + -,t) dxdkdadt 
Jr V2fc e 
,t ,2tt , ^ ^x,^,^,*) sine* i 

/ / / — : -= {o a ip) (x, fc, a H — ,t) dx.dk da dt 

Jo Jo Jr \2k e 



/o Jo 

r T /.2ir 



x (^) (* + 721 sin a, * - V2* cos a , « + J, t ) .x dfc da * 

'-" ' sma 
£ e ,i(x,t) 



o Jo Jr v2fc 
x (^) (* + VZk sin a, * 2 - V2* cos a, *, a + J, t ) cix dfc cia * , 



and 

cT r 2-n 



n( <9 2 c/> e (x, k,a , t) ip(x., k, a, t) dxdk da dt 
Jr e 

= f f f V2k E et i(x,k,a,t) cosa (d a -)p)(x,k,a + -,t) dxdk dadt 
Jo Jo Jr e 



to Jo Jr 

+ [ [ [ V / 2fc-B e ,2(x, k, a, t) sin a (d a ip)(x, k, a + -, t) dxdk da dt 
Jo Jo Jr e 



lo Jo Jr 



x {d a tj))(x\ + V2k sin a, x 2 — V2k cos a, k, a + -, t) di.dk da dt 

e 

rT r 2ir 



/ E ei2 {x,t)V2k sin a 
Jr 



/o Jo Jr 

x (d a ip)(x-L + V2k sin a, x 2 — cos a, fc, a + -, t) dxdk da dt , 
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(3.47) 



/o Vo Jr v2fc 



S e ,i(x,i)\/2fc cos a ( 3 - 49 ) 



so we have 

pT r 2ir 



lim f [ [ d\ a (j) t (x,k,a ,t) tp(x, k,a,t) dxdk da dt 

: ^°Jo Jo Jr ' 6 



T /.2-7T 

' - ' ' - , x cos a 
E 2 (x,t) 



lo Jo Jr v2fc 

f-2-K 

X 



r-ZTt 

/ {d a ip){xi + V2k sin a, x 2 — V2fc cos a, k,a + r, t) drdxdk dadt 
Jo 

[ r f E^t) 

Jo Jo Jr 



'2k 

r 2ir 



r — 

x / (d a tp)(xi + V2k sin a, £2 — V2k cos a, k, a + r, t) dr dii.dk da dt 
Jo (3.50) 

/ E (~ t) C ° Sa 
Jr V2k 



r-2-K 

x / (d T ip)(xi + V2k sina,X2 — V2k cos a, k,a + r, t) drd5s.dk dadt 
Jo 



T r 2tr 



~ , sin a 



lo Jo Jr y2k 

I>2tt 

x / (d T ip){xi + V2fc sina,X2 — V2k cos a, k,a + t, t) dr dxdk dadt 
Jo 

= 0, 



and 



T p 2-k 



lim III dl a e (x, k, a ,t) tp(x, k,a,t) dxdk da dt 

e ^°Jo Jo Jr ' e 



T r 2-K r 

J Ei(x,t) V2k cos a 



/o Jo Jr 

v2tx 

< ' 

io 

r T i-2ix 



/ {d a ^){xi + V2k sina,x 2 — V2k cos a, fc,a + r,t) dr dxdk dadt 
Jo 

+ [ [ [ E 2 (x,t) V2k sin a 

Jo Jo Jr 

r 27T _ _ 

x / (d a xp)(xi + V2k sina,X2 — V2k cosa,k,a + r,t) dr dxdk dadt 
Jo (3.51) 

/ Si(x,t) V2fc cos a 
Jr 

/• 27r _ _ 

x / (d T ip)(xi + V2k sin a, 2; 2 — V2k cos a, fc,a + r, i) dr dxdk dadt 
Jo 

/ £ 2 (x,t) \/2fc sin a 
Jr 

r 27T 

/ (d T ip)(xi + V2k sina,X2 — ^2k cos a, k,a + r,t) dr dxdk dadt 
Jo 



r27T 

- 

0. 



Then we deduce that 



V (fe , Q )9 Q ^ ${x,k,a-T,T,t)d^J =0, (3.52) 
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( r 27r \ ■ 

and that the function (x, k,a) i— > d a \ / <£>(x, fc, a — r, r, t) dr is constant on f2. Then, 

for any (x, k, t), the function an-/ <I>(x, fc,a — T,T,t)dr is monotonic and periodic on 

Jo 

[0, 27r] and then must be constant. As a consequence, we finally obtain 

<9 Q ^ <P(x,k,a-T,T,t)dT^j = 0. (3.53) 



□ 



Consequently, if we consider the hypothesis (|3.41[) alongwith those from Theorem [31 we 
can write 



/ <£>(x, k, a — r, r, t) dr = / <£>(x, k, — r, r, i) dr . 
Jo Jo 

Then, the two-scale limit model Q3.13P can be reduced to 
( d t G + (£ 2 ) d Xl G— {£{) d X2 G+ (T a ) d a G = , 

G(x, fc, a, 0) = —f°(xi — \/2k sin a, x 2 + \/2k cos a, \/2fc cos a, -\/2fc sin a) , 
-Vxl(x,T,t) = £(x,T,t), 

— Ai$(x,r, t)= / G(xi + \/2fc sina,X2 — V2fc cos a, k,a + r, i) dkda 
Js 

-^n e (x), 

, , i?2(x, fc, a, r, £) cosa — £i(x, fc, a, r, i) sina 

a( x i k, T, t) = 



(3.54) 



(3.55) 



V2fc 

and the notation (•) stands for 

(u)(x.,k,t)= / u(x, k, —T,r,t) dr . (3.56) 
./o 

Theorem 7. Under the hypotheses of Theorem^ and \3.4-l\j , there exists a function g defined 
on r x [0, T] such that 

f(x,k,a,t) = ^-g(x,k,t), V(x,fc,a,i) Gfix [0,T], (3.57) 
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and verifying 

d t g + (£2) d Xl g - (Si) d X2 g = 

,-2n 



/>Z7T 

g(x, k, 0) = / f°(xi — V2k sina,X2 + V2k cos a, y2k cos a, V2k sin a) da , 
J 

(S)hc,k,t) = — / E(xi - \/2fc sina,a; 2 + \p2k cosa,t)da, (3.58) 
27r Jo 

-V^(x,i)=E(x,t), 

— Aj^(x, i) = — / 5(^1 + v^fc sin a, X2 — V%k cos a, k, t) dk da — n e (x) . 
2tt J 5 

Proof of Theorem [?| Since F and G are linked by the relation p. lip , / does not depend on 
a. Indeed, we have 

d a f(x,k 7 a 7 t)~d a / F(x,k 7 a 7 T,t) dr 
Jo 

p2ir 

— d a / G(x,k,a + T,t)dr 
Jo 

p2lT 

= d a G(x, k,r,t) dr 
Jo 



(3.59) 



= 0. 

Then, if we integrate (|3.55l a) in a, we obtain 

d t f+(S 2 )d x J-(S 1 )d x J = 0, (3.60) 

and the equation (|3.58l a). 

By integrating the initial condition (|3.55l b) in a and dividing it by 2-7T, we obtain the 
initial condition (13.581 ^. 

Since (S ) does not depend on a, we can integrate it in a and divide it by 2tt: then, we 
have 

^ />2"7T />27T 

(S)(x,k,t) = — / / £(x, fc,a — T,r,t) dr da 
27r Jo Jo 

= — I / £(x,k 7 a 7 T 7 t)dTda (3.61) 
2tt Jo Jo 

1 />2tt />27r 

= — / / ^(xi ~ \ // 2k sma,X2 + V2k cos a,t) da . 
27r Jo Jo 

By using similar techniques, we integrate (|3.55l d) in r and we obtain 

— Ax0(x, t) = — [ g( xi +V2k sina,X2 — V^k cos a, k,t)dk da — n e (x) , (3.62) 
2tt J 5 

Finally, we integrate p. 551 c) in r and we obtain (|3.58l d). which concludes the proof. □ 

We can remark that, even if they are based on the non-physical hypothesis p.4ip for 
(E e ) e> o, Bostan's results on the mathematical justification of the 2D finite Larmor radius 
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approximation have been improved: indeed, we have proved that the hypotheses of Theorem 
[3] coupled with (|3.41|) are sufficient to prove Theorem [71 Furthermore, this result has been 
generalized to a non-periodic initial distribution /o and to a non-uniform electron density 

4 Conclusions and perspectives 

In a first part, we recalled a two-scale convergence result for a 2D Vlasov-Poisson model for 
a charged particle beam which is due to Frenod and Sonnendriicker. After adapting it to 
plasma modeling by adding an electron density h e and some reasonable compatibility con- 
ditions, we recalled that such a result trivially implies, up to some subsequences, the weak-* 
convergence of (/ e ,E e ) to a couple (/,E). Then we compared this conclusion to Bostan's 
results which was established in [5] under stronger assumptions. 

In a second part, we introduced a new set of variables involving the guiding-center po- 
sition coordinates and the transverse part of the kinetic energy which is denoted with k. 
After rewriting the Vlasov equation in these new coordinates, we established a two-scale 
convergence result by using Frenod & Sonnendriicker's assumptions and the compatibility 
conditions for n e which were added in Theorem [TJ Then, we proved that (/ e ,E e ) weakly-* 
converges to a couple (/, E) and, under an additional non-physical property for (E e ) e> o, 
we established a system of constraints satisfied by (/, E) through a few computations lines, 
remarking that this limit model is exactly the 2D finite Larmor radius model which was 
studied in [5]. 

The first main remark we can do about the present work concerns the formulation of the 
two-scale limit model in canonical gyrokinetic coordinates under hypothesis (|3.4ip : indeed, 
the transport equations (|3.55l a) and (I3.58l a) indicate that k is an adiabatic invariant not 
only for the weak-* limit system but also for the two-scale limit system. Since k is the dimen- 
sionless transverse kinetic energy, this remark confirms previous results about the adiabatic 
invariant property of the magnetic momentum (see Littlejohn [21], Lee [HI [20], Brizard el 
al. [6] [7j, Grandgirard et al. [IT]). 

The second remark which can be done is about the mathematical results which have 
been established in this paper: by linking previous two-scale and weak-* convergence under 
common assumptions, i.e. the assumptions from Theorem [3] coupled with f|3 .41 [) . we have 
improved Bostan's results by generalizing them to some cases involving a non-uniform elec- 
tron density or a non-periodic initial distribution /o. 

Since the canonical gyrokinetic variables allow us to simplify the formulation of the two- 
scale limit model and, under (|3.41l) . allow us to highlight the adiabatic invariant property of 
the magnetic momentum, they may be a useful tool for a mathematical justification of the 
full 3D finite Larmor radius model. 

From a numerical point of view, the model (|3.13|) can be used to build a two-scale 
numerical method in order to simulate the high frequency oscillations of the solution (/ e , E e ) 
of the Vlasov-Poisson model (12 . 16[) . such as it has been done before for low Mach number 
problems (see [9]), charged particle beams problems (see [IT] and [24]) or drift problems in 
the ocean (see [2]). As an example, a numerical method based on the computation of the 
characteristics associated with the limit transport equation will be simpler to be developed 
on the formulation (I3.13l a) than the formulation (|2.21l a): indeed, since we have 

d Xl {£ 2 ) - d X2 (£i) = d k {F k ) + d a {F a ) = , (4.63) 

the transport equation (|3.13l a) can support a time splitting without losing its conservation 
property (see [4j). 



19 



Acknowledgments: Wc would like to thank Daniel Han-Kwan for a helpful remark he 
did and which allowed us to correct a mistake. 



References 

[I] R. Adams, J. Fournier, Sobolev Spaces, Academic Press (1975 - 2 nd ed. in 2003). 

[2] P. Ailliot, E. Frenod, V. Monbet, Long term object drift in the ocean with tide and 
wind, Multiscale Model. Simul. 5-2 (2006), 514-531. 

[3] G. Allaire, Homogenization and two-scale convergence, SIAM J. Math. Anal. 23-6 
(1992), 1482-1518. 

[4] N. Besse, E. Sonnendrucker, Semi-Lagrangian schemes for the Vlasov equation on 
an unstructured mesh of phase space, J. Comp. Phys. 191 (2003), 341-376. 

[5] M. Bostan, The Vlasov-Poisson system with strong external magnetic field. Finite Lar- 
mor radius regime, Asymptot. Anal. 61-2 (2007), 91-123. 

[6] A. Brizard, Nonlinear gyrokinetic tokamak physics, PhD thesis of Princeton University 
(1990). 

[7] A. Brizard, T.-S. Hahm, Foundations of nonlinear gyrokinetic theory, Rev. Mod. Phys. 
79 (2007), 421-468. 

[8] D.-H. Dubin, J. -A. Krommes, C. Oberman, W.-W. Lee, Nonlinear gyrokinetic 
equations, Phys. Fluids 26-12 (1983), 3524-3535. 

[9] E. Frenod, A. Mouton, E. Sonnendrucker, Two-scale numerical simulation of the 
weakly compressible ID isentropic Euler equations, Numer. Math. 108-2 (2007), 263-293. 

[10] E. Frenod, P. -A. Raviart, E. Sonnendrucker, Two-scale expansion of a singularly 
perturbed convection equation, J. Math. Pures Appl. 80-8 (2001), 815-843. 

[II] E. Frenod, F. Salvarani, E. Sonnendrucker, Long time simulation of a beam in 
a periodic focusing channel via a two-scale PIC-method, Math. Models Methods Appl. 
Sci. 19-2 (2009), 175-197. 

[12] E. Frenod, E. Sonnendrucker, Homogenization of the Vlasov equation and of the 
Vlasov-Poisson system with a strong external magnetic field, Asymptot. Anal. 18-3-4 
(1998), 193-214. 

[13] E. Frenod, E. Sonnendrucker, Long time behavior of the Vlasov equation with a 
strong external magnetic field, Math. Models Methods Appl. Sci. 10-4 (2000), 539-553. 

[14] E. Frenod, E. Sonnendrucker, The finite Larmor radius approximation, SIAM J. 
Math. Anal. 32-6 (2001), 1227-1247. 

[15] F. Golse, L. Saint-Raymond, The Vlasov-Poisson system with strong magnetic field, 
J. Math. Pures Appl. 78 (1999), 791-817. 

[16] F. Golse, L. Saint-Raymond, The Vlasov-Poisson system with strong magnetic field 
in quasineutral regime, Math. Models Methods Appl. Sci. 13-5 (2003), 661-714. 

[17] V. Grandgirard, M. Brunetti, P. Bertrand, N. Besse, X. Garbet, P. Ghen- 
drih, G. Manfredi, Y. Sarazin, O. Sauter, E. Sonnendrucker, J. Vaclavik, 
L. VlLLARD, A drift-kinetic semi-lagrangian J^D code for ion turbulence simulation, J. 
Comp. Phys. 217 (2006), 395-423. 

[18] D. Han-Kwan, The three-dimensional finite Larmor radius approximation, Asymptot. 
Anal. 66-1 (2010), 9-33. 



20 



[19] W.-W. Lee, Gyrokinetic approach in particle simulation, Phys. Fluids 26-2 (1983), 
555-562. 

[20] W.-W. Lee, Gyrokinetic particle simulation model, J. Comp. Phys. 72-1 (1987), 243- 
269. 

[21] R.-G. Littlejohn, A guiding center Hamiltonian : A new approach, J. Math. Phys. 
20-12 (1979), 2445-2458. 

[22] J.-L. Lions, Quelques methodes de resolution de problemes aux limites non lineaires, 
Dunod, Gauthicr-Villars (1969). 

[23] A. Mouton, Approximation multi-echelles de V equation de Vlasov, these de l'Universite 
de Strasbourg, Ed. Universitaires Europeennes, TEL-00411964 (2009). 

[24] A. MOUTON, Two-scale semi-lagrangian simulation of a charged particle beam in a 
periodic focusing channel, Kinet. Relat. Models 2-2 (2009), 251-274. 

[25] G. Nguetseng, A general convergence result for a functional related to the theory of 
homogenization, SIAM J. Math. Anal. 20-3 (1989), 608-623. 



21 



